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Abstract 

We consider statistical models driven by Gaussian and non-Gaussian self-similar pro- 
cesses with long memory and we construct maximum likelihood estimators (MLE) for the 
drift parameter. Our approach is based in the non-Gaussian case on the approximation by 
random walks of the driving noise. We study the asymptotic behavior of the estimators 
and we give some numerical simulations to illustrate our results. 
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1 Introduction 

The self-similarity property for a stochastic process means that scaling of time is equivalent to 
an appropriate scaling of space. That is, a process (Yt)t>o is selfsimilar of order H > if for 
all c > the processes (Y c t)t>o and (c H Yt)t>o have the same finite dimensional distributions. 
This property is crucial in applications such as network traffic analysis, mathematical finance, 
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astrophysics, hydrology or image processing. We refer to the monographs [T], [I] or [TU] for 
complete expositions on theoretical and practical aspects of self-similar stochastic processes. 

The most popular self-similar process is the fractional Brownian motion (fBm). Its 
practical applications are notorious. This process is defined as a centered Gaussian process 
(BP)t>o with covariance function 

R H (t, s) := E(B?B?) = ~ (t 2H + s 2H -\t- s\ 2H ) , t,s>0. 

It can be also defined as the only Gaussian self-similar process with stationary increments. 
Recently, this stochastic process has been widely studied from the stochastic calculus point of 
view as well as from the statistical analysis point of view. Various types of stochastic integrals 
with respect to it have been introduced and several types of stochastic differential equations 
driven by fBm have been considered (see e.g. [8], Section 5). Another example of a self- 
similar process still with long memory (but non-Gaussian) is the so-called Rosenblatt process 
which appears as limit in limit theorems for stationary sequences with a certain correlation 
function (see [2], p3]). Although it received a less important attention than the fractional 
Brownian motion, this process is still of interest in practical applications because of its self- 
similarity, stationarity of increments and long-range dependence. Actually the numerous uses 
of the fractional Brownian motion in practice (hydrology, telecommunications) are due to these 
properties; one prefers in general fBm before other processes because it is a Gaussian process 
and the calculus for it is easier; but in concrete situations when the Gaussian hypothesis is not 
plausible for the model, the Rosenblatt process may be an interesting alternative model. We 
mention also the work |14j for examples of the utilisation of non-Gaussian self-similar processes 
in practice. 

The stochastic analysis of the fractional Brownian motion naturally led to the statistical 
inference for diffusion processes with fBm as driving noise. We will study in this paper the 
problem of the estimation of the drift parameter. Assume that we have the model 

dX t = 9b(X t )dt + dBf, t G [0, T] 

where (-Bf0ie[0,T] i s a fractional Brownian motion with Hurst index H 6 (0, 1), b is a deter- 
ministic function satisfying some regularity conditions and the parameter 9 £ R has to be 
estimated. Such questions have been recently treated in several papers (see [6], [16] or [12]): in 
general the techniques used to construct maximum likelihood estimators (MLE) for the drift 
parameter 9 are based on Girsanov transforms for fractional Brownian motion and depend on 
the properties of the deterministic fractional operators related to the fBm. Generally speaking, 
the authors of these papers assume that the whole trajectory of the process is continuously 
observed. Another possibility is to use Euler-type approximations for the solution of the above 
equation and to construct a MLE estimator based on the density of the observations given "the 
past" (as in e.g. |9j], Section 3.4, for the case of stochastic equations driven by the Brownian 
motion). 

In this work our purpose is to make a first step in the direction of statistical inference 
for diffusion processes with self-similar, long memory and non-Gaussian driving noise. As 
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far as we know, there are not many result on statistical inference for stochastic differential 
equations driven by non-Gaussian processes which in addition are not semimartingales. The 
basic example of a such process is the Rosenblatt process. We consider here the simple model 



where (^t )te[o,T] is a Rosenblatt process with known self-similarity index H £ (2, 1) (see 
Sections [3] and Appendix for the definition) and a £ M. is the parameter to be estimated. We 
mention that, since this process is not a semimartingale, it is not Gaussian and its density 
function is not explicitly known, the techniques considered in the Gaussian case cannot be 
applied here. We therefore use a different approach: we consider an approximated model in 
which we replace the noise Z H by a two-dimensional disturbed random walk Z H,n that, from 
a result is [15], converges weakly in the Skorohod topology to Z H &s n — > 00. Note that this 
approximated model still keeps the main properties of the original model since the noise is 
asymptotically self-similar and it exhibits long range dependence. We then construct a MLE 
estimator (called sometimes in the literature, see e.g. [9] "pseudo-MLE estimator") using an 
Euler scheme method and we prove that this estimator is consistent. Although we have not 
martingales in the model, this construction involving random walks allows to use martingale 
arguments to obtain the asymptotic behavior of the estimators. Of course, this does not solve 
the problem of estimating a in the standard model defined above but we think that our approach 
represents a step into the direction of developing models driven by non-semimartingales and 
non-Gaussian noises. 

Our paper is organized as follows. In Section[2]we recall some facts on the pseudo MLE 
estimators for the drift parameter in models driven by the standard Wiener process and by 
the fBm. We construct, in each model, estimators for the drift parameter and we prove their 
strong consistency (in the almost sure sense) or their L 2 consistency under the condition a > 1 
where N a is the number of observations at our disposal and the step of the Euler scheme 
is jf. This condition extends the usual hypothesis in the standard Wiener case (see [2j see 
also [9], paragraph 3.4). Section [3] is devoted to the study of the situation when the noise is 
the approximated Rosenblatt process; we construct again the estimator through an inductive 
method and we study its asymptotic behavior. The strong consistency is obtained under similar 
assumptions as in the Gaussian case. Section 0] contains some numerical simulations and in the 
Appendix we recall the stochastic integral representations for the fBm and for the Rosenblatt 
process. 

2 Preliminaries 

Let us start by recalling some known facts on maximum likelihood estimation in simple stan- 
dard cases. Let (Wt)tg[o,T] be a Wiener process on a classical Wiener space (fi,^", P) and let 
us consider the following simple model 



X t = at + Z i 



't 1 





(1) 
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with T > and assume that the parameter a G M has to be estimated. One can for example 
use the Euler type discretization of dH) 



:= Y tj+1 = Y tj + aAt + W tj+1 -W tj , j = 0,...,N-l, 

with Yi =Yq = and At = tj+i — tj the step size of the partition. In the following, we denote 
Yf. =Yj. In the following fz denotes the density of the random variable Z. 

The conditional density of with respect to Yi,...Yj is, by the Markov property, 
the same as the conditional density of with respect to Yj and since Wt j+1 — has the 
normal law N(0, At), this density can be expressed by 

1 / 1 (y J+1 - yj - aAtf 



/ W ^(W^) = 7 =exp( 



y/2ir(At) V 2 At 
We easily obtain the likelihood function of the observations Y\, . . . , Yat 



JV-l „ / „ N-l 

L(a,yi, ...,y N )= f Yl (yi) fY j+1 /Y s (yj+i/Vj) 

j=i 



{2*{At)) N / 2 P l 2^L At 



and this gives a maximum likelihood estimation of the form 

N-l 

NAt 



3=0 

and the difference ajy — a can be written as 

N-l 



i=o 

Then 

,2 1 



E ajv ~~ a \ 



NAt' 

and this converges to zero (that is, the estimator is Inconsistent) if and only if 

NAt ->■ oo, N -> oo. (2) 
Note that the partition = with j = 0, . . . N does not satisfy ([2]). 

Remark 1 Under condition |||) we get, fry t/ie strong law of large numbers, the almost sure 
convergence of the estimator a at to the parameter a. 
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Remark 2 We need in conclusion to consider an interval between observation of the order 
= -j^ with a < 1 to have Equivalently, if we dispose on N a observations with a > 1, 
i.e. T > N"" 1 , and the interval At is of order jj, condition (djj still holds. Using this fact, 
we will denote in the sequel by N a the number of observations and we will use discretization 
of order -k of the model. 



Next, let us take a look to the situation when the Brownian motion W is replaced by 
a fractional Brownian motion B H with Hurst parameter H £ (0, 1). As far as we know, this 
situation has not been considered in the literature. The references [6] or [16] uses a different 
approach, based on the observation of the whole trajectory of the process Yj below. The model 
is now 

Y t = at + Bf, te[0,T\, (3) 

and as before we aim to estimate the drift parameter a by assuming that H is known and on 
the basis on discrete observations Yj, . . . , Y/v<* (the condition on a will be clarified later). We 
use the same Euler type method with tj = and we denote Y tj = Yj. 

We can easily found the following expression for the observations Yj, j = 1, . . . , N a , 

We need to compute the density of the vector (Yj , . . . , Yj^a ) . Since the covariance matrix of 
this vector is given by V = (Ti j)i t j=i n<* with 

Tij = Cov (BfBf 

\ N N 

the density of (Yj, . . . , Yn<*) will be given by 

(2vr)-^ - T L= exp ( -- U , y N « - N a -Y I?" 1 U - -, . . . , y N « - N a - 

VdrtT V 2 ^ N N ' v N' ,y N 

and by maximizing the above expression with respect to the variable a we obtain the following 
MLE estimator 

a N - N Na : (4) 



and then 

\^N a --p-1 r>H 
Z^i,j=l J i,j D ±_ 
a N - a = N a N , (5) 

where the Ty} are the coordinates of the matrix 
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Thus 



E\a N -a\ 2 = N 



, '■>' >J t \ N N / 



Note that 

E ^r^(5|B|)= e ^fe* = E i^=E^- 

i,j,k,l=l j,k,l=l \i=l / j,k,l=l i>'=l 

and consequently 

it. I - i2 _ Ar 2 1 _ A7-2-2H 1 

aiv — a — iv —jga r = 



EAfa . v~-7V Q • ■ -1 



where ™^ are the coefficients of the matrix M 1 with M = (571^)^=1,. ..,jv-, w-ij = \{i 2H + 
j 2iT — |i — j| 2H ). Let x be the vector (1, 2, ... , iV Q ) in R^". We use the inequality 

x t M~ 1 x > ^2 
A 

where A is the largest eigenvalue of the matrix M. Then we will have 

E|aAr — a| < iv 2 . 

\\ x II2 

Since l 2 + 2 2 + . . . p 2 = p ( p+1 ^ 2p+1 ) we see that ||x|| 2 behaves as iV 3 " and on the other hand 
(by the Gersghorin Circle Theorem (see [5], Theorem 8.1.3, pag. 395) 

N a 



A < max V \m u \ < CN a{ - 2H+l \ 

~ i=l,... JV- 1 ' ~ 



1=1 

with C a positive constant. Finally, 

E \a N - a\ 2 < CN 2 - 2H - 3a+a ( 2H -V = cw(2-2fi)(--«) ( 6 ) 

and this goes to zero if and only if a > 1. 

Let us summarize the above discussion. 

Proposition 1 Let (-B^)t-[o,T] be a fractional Brownian motion with Hurst parameter H 6 
(0, 1) and let a > 1. T/ien the estimator ^ is LP consistent for any p > 1. 
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Proof: Since for every n the random variable a n — a is a centered random variable, it holds 
that, for some positive constant depending on p 

p 

E \a N - a\ p < c p (E \a N - a| 2 ) 2 < ^JVJ»(i-#)(i-«) 
and this converges to zero as zero since a > 1. ■ 

It is also possible to obtain the almost sure convergence of the estimator to the the 
true parameter from the the estimate ©. 

Proposition 2 Let {B^) t ^ T ^ be a fractional Brownian motion with Hurst parameter H € 
(0, 1) and let a > 1. Then the estimator ^ is strongly consistent, that is, a n — a -^ n -+oo P 
almost surely. 

Proof: Using Chebyshev's inequality, 

P (\a n -a\> N-P} < CpN^N^ 1 - 11 ^ 1 -^. 

In order to apply the Borel-Cantelli lemma, we need to find a strictly positive (5 such that 
J2n>i N^N^ 1 - 11 ^ 1 -^ < oo. One needs pf3 + (1 - H){1 - a)p < -1 and this is possible if 
and only if a > 1. ■ 

Remark 3 The fact that the restriction a > 1 is interesting because it justifies looking for 
another type of estimator for which no such restriction holds. But for the approach using the 
Euler scheme discretization, this restriction is somehow expected because it appear also also in 
the Wiener case (see J2J/j. 

Remark 4 Let us also comment on the problem of estimation of the diffusion parameter in 
the model $3$. Assume that the fractional Brownian motion B H is replaced by o~B H in $3$, 
with a € R. In this case it is known that the sequence 

AT-l 

N 2H-iy ( Yi+1 - Yi ) 

^ \ — nJ 

i=0 

converges (in L 2 and almost surely ) to a 2 . Thus we easily obtain an estimator for the diffusion 
parameter by using such quadratic variations. The above sequence has the same behavior if we 
replace the fBm B H by the Rosenblatt process Z H because the Rosenblatt process is also self- 
similar with stationary increments and it still satisfies E \ Zp — = \t — s\ 2H for every 
s,t G [0, T] (see Section 4 for details). For this reason we assume throughout this paper that 
the diffusion coefficient is equal to 1. 
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3 MLE and random walks in the non- Gaussian case 



We study in this section a non-Gaussian long-memory model. The driving process is now 
a Rosenblatt process with selfsimilarity order H E This process appears as a limit 

in the so-called Non Central Limit Theorem (see [2] or |13|). It can be defined through its 
representation as double iterated integral with respect to a standard Wiener process given by 
equation (|13p in the Appendix. Among its main properties, we recall 

• it exhibits long-range dependence (the covariance function decays at a power function at 

zero) 

• it is iJ-selfsimilar in the sense that for any c> 0, (Z H (ct)) =^ (c H Z H (t)), where " = (d) " 

means equivalence of all finite dimensional distributions ; moreover, it has stationary 
increments, that is, the joint distribution of (Z H (t + h) — Z (h),t £ [0,T]) is independent 
of h > 0. 

• the covariance function is 

E(Zf Zf ) = \ {t 2H + s 2H -\t- s\ 2H ) , s, t € [0, T] 

and consequently, for every s,t E [0, T] 

E\Z?-Z?\ 2 = \t-s\ 2H 

• the Rosenblatt process is Holder continuous of order 5 < H 

• it is not a Gaussian process; in fact, it can be written as a double stochastic integral of a 

two-variable deterministic function with respect to the Wiener process. 

Assume that we want to construct a MLE estimator for the drift parameter a in the 

model 

Y t = at + Z t H , te[0,T\. 

We first note that the Rosenblatt process has only been defined for the self-similarity order 
H > i and consequently we will have in the sequel always H E (1, 1). The approaches used 
previously do not work anymore because the Rosenblatt process is still not a semimartingale 
and moreover, in contrast to the fBm model, its density function is not explicitly known 
anymore. The method based on random walks approximation offers a solution to the problem 
of estimating the drift parameter a. We will use this direction; that is, we will replace the 
process Zh by its associated random walk 

zr= E ^ 2 /* / # f(M«,,W^M^ te[o,T] (7) 
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where the £j are i.i.d. variables of mean zero and variance one and the deterministic kernel F 
is defined in Appendix by (I14|) , It has been proved in |T5] that the random walk (|7|) converges 
weakly in the Skorohod topology to the Rosenblatt process 



We consider the following discretization of the Rosenblatt process 

(Z^ H ), j = 0,...N a , 

N 

where for j / 1 Z / is given by (J7|) and we set Z 1 ^ = £i/N H . With this slight modification, 

N AT 

the process (Z^) still converges weakly in the Skorohod topology to the Rosenblatt process 

JV 

Z H . We will assume as above that the variables £j follows a standard normal law iV(0, 1). 

Concretely, we want to estimate the drift parameter a on the basis of the observations 

Y tj+1 = Y h + a(t j+1 - tj) + {Z*£ - zf: N 

where tj = -b, j = 0, . . . , N a and Yq = 0. We will assume again that we have at our disposal 
a number N a of observations and we use a discretization of order of the model. Denoting 
Yj := Y tj , we can write 

Y j+1 = Y j + ^ + (z^-Z^, j = 0,...,N*-l. 

Now, we have 

z £f - Z ± N = fj(£h ■ • • >£j) + 9j(tl, 

N N 

where f =0,h = = -£i/N H and for j > 2 

U ti) = N E (f (^T' u,v^-F (jj,u, v^j ) dvduj 

and g = 1/N H for j > 1 

9j =#(&,..., go = 2Ar E ^ F n ' w ) dudu J 

Finally we have the model 

Yj+^Yj + ^ + fj+gjSj+v (8) 
In the following, we assume that . . . , £ n , . . .) £ where 

s = n i >i{ 5j (a,---,ej)^o}. 

The event £ satisfy ¥(B) = 1. 
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Remark 5 Note that on the event B, conditioning with respect to the - • • , £j is the same 
as conditioning with respect to the Y± , ■ ■ ■ ,Yj. In fact, we have 

Y 1 = a/N + ^/N H 

Yj = Y j - 1 +a/N + f- 1 (Z 1 ,...,£ j - 1 )+£ j g j - 1 (Z 1 ,...,Z j -i), j>2. 

Since on B, for all j > 2, <?j_i(£i, . . . , / 0, i/te too a-algebra <r(£i, . . . , cmd 

£7(li, . . . sato/t/ 

. . . , nB = o-(F!, . . . , y,_i) n b. 



Then, given the random variable Yj+i is conditionally Gaussian and the 

conditional density of Yj + \ given Y±, . . . ,Yj can be written as 

fY j+1 /Y 1 ,...,Y j {yj+i/yi,---,yj) = -^=? ex P ^ — • 



j j 

The likelihood function can be expressed as 



:i 



L(a,yi,...,y Na ) = f Yl (yi)fY 2 /Y 1 {y2/yi) • • • fY Na /Y 1 ,...,Y Na _ 1 (VN<*/Vi, ■■■,VN°-i) 



JVa - 1 1 / 1 ..._„/JVT_ t.\2 



l(Vj+i-yj-a/N - fj) 



o a/27T3 



- ex P , o 

2 V 2 gr 



2 



By standard calculations, we will obtain 



T 

^3 

and since Yj+i — Yj — fj = + gj^j+i we obtain 



& N = ^N--l 1 3 ( 9 ) 



Let us comment on the above expression. Firstly, note that since z2j=o ~g~? ^ s no ^ 
deterministic anymore, the square mean of the difference ajy — a cannot be directly computed. 
Secondly, if we denote again by 

M-\ 



3=0 
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this discrete process still satisfies 

E(A m+ i/Fm)=A m , VM>1 

where Tm is the cr-algebra generated by £1, . . . ,£m- However we cannot speak about square 
integrable martingales brackets because A 2 M is not integrable (recall that the bracket is defined 

in general for square integrable martingales). In fact, E(^4| f ) = Ejf=o^ fs^J an< ^ * S 

finite in general because gj is a normal random variable. We also mention that, in contrast to 
the Gaussian case, the expectation of the estimator is not easily calculable anymore to decide if 
([9]) is unbiased. On the other hand, from the numerical simulation it seems that the estimator 
is biased. 

Nevertheless, martingale type methods can be employed to deal with the estimator ([9]). 
We use again the notation 

M-l 
= E 72- 

j=o y i 

The following lemma is crucial. 
Lemma 1 Assume that a > 2 — 2H and let us denote by 

N a -l 

i i 

Tn := N2 1^ 72- 



Then T/v N ~*°°) oo almost surely. Denote 



N 2 — g- 

3=0 3 3 



s 1-7 

N - 



Then there exists < 7 < 1, such that Un — ~*°° > 00 almost surely. 

Proof: Let prove the convergence for Tjy. We will use a Borel-Cantelli argument. To this 
end, we will show that 

P (T N < iV 5 ) < 00 (10) 

N>1 

for some 5 > 0. 

Fix < 5 < a - (2 - 2H). We have 



P (t n < N 5 ^ 




N a — 1 1 
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We now choose a p integer large enough such that | < a — (2 — 2H) — 5 and we apply the 
Markov inequality. We can bound P (T N < N 5 ) by 

P (t n < N 5 ^ < N^ 2+S ^E 



2^7=0 a? 



and using the inequality "harmonic mean is less than the arithmetic mean" we obtain 



(T N < iV 5 ) < N^ p N- 2ap K Yj 9j 



j=0 

Note that (the first inequality below can be argued similarly as in 

2 



1 



N 2H 



> E 



iyH,N ryH^ 

JV N 



it holds that 



Hal) :=c ,<]v^- 

Since £j = y/cjXj where X,- ~ N(0, 1), by © 



(N a -1 
i=o 

< N (2+5-2a-2H)p N ap _ N p(2+8-a-2H) 

and thus relation (fT0|) is valid. 

The convergence for Un can be obtained in a similar way with < 7 < 1 such that 
a(2-j)-2+2Hj > 0, < 5 < a(2-j)-2+2Hj andp such that 1/p < a(2-~f)-(2-2Ha)-5. 
U 

We use the notation 

and 

(A) Af+ i - (A) M 
Bm := — — • 

\ A /M+1 

Note that Vm and 5j| are Tm adapted. 
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We recall the Robbins-Siegmund criterium for the almost sure convergence which will 
play an important role in the sequel (see e.g. [3j, page 18): let (Vn)n, (Bn)n be J^n adapted, 
positive sequences such that 

E(V N +i/T N ) <V n + B n a.s. . 

Then the sequence of random variables (Vn)n converges as N — > oo to a random variable Vx, 
almost surely on the set {J2n>i < 00 }• 

Lemma 2 The sequence Vm converges to a random variable almost surely when M — > oo. 
Proof: We make use of the Robbins-Siegmund criterium. It holds that 

( a 2 . 

E(V m +x/Fm)=® 




1 

(A) 

M+l 

' (A) 

M+l 

< Vm + %■ 



< [a m+1 /t m ) = —T^r {A* + ( a )m+i - (A) M ) 



By Lemma [J the sequence {A)m converges to oo and therefore 

/•oo 

y"B Na <C+j x- 1 -~ 1 ds<oo. 

N Jl 

Once can conclude by applying the Robbins-Siegmund criterium. 

We state now the main result of this section. 
Theorem 1 The estimator is strongly consistent. 
Proof: We have 

(a N - a) 



2 V/v« Vn° 



and we conclude by Lemmas Q] and [2j ■ 

Comment: The L 1 (or LP) consistence of the estimator ([9]) is on open problem. Note that 
its expression is a fraction whose numerator and denominator are non-integrable involving 
inverses of Gaussian random variables. The natural approaches to deal with do not work. 
A first basic idea when is to try to apply the Holder inequality for the product FG with 
F = y^^_ n _1 ^ii- and G = y^^_ n _1 \- But when we do this, we are confronted with the 



moments of the random variable \F\ and it has not moments. Indeed, the second moment 
is E^feo" 1 "Tf w hich is obviously infinite. Even E|jP| is infinite. Another way is to use the 
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inequality (X^ILi l°A|) 2 — (^27=1 l a «| 2 ) {Yli=i l^«| 2 )- l n this wa y we av °id the non-integrability 
of the random variable F defined above. We get 



N a -1 




and then 



\a-N — a\ < N 



N a -1 1 

i=° 5= 



and by saying that the harmonic mean is less that the arithmetic mean (we have the feeling 
that it is rather optimal in this case) 



(yv Q — 1 \ 2 I N a -i 
j=0 J \ j=0 

We tried now to apply Holder for E|oat — a\ and this gives 



l-a 




(11) 



E|aAr - a\<N 



But now S^Lo 1 is a chi-square variable with iV a degree of freedom, its p moment behaves 

,Ara_1 " 2 i 2 i is of order less than N- a 2N~ H 



as iV Qp and it can be also proved that I E ( X^/=o 5 



the paper it is proved in the case q = 2) which unfortunately gives X^o"* — ^ 
One needs to compute directly the norm of \cin — o>\ in (|lip by using the distribution of the 
vector ((71, ... , But it turns out that this joint distribution is complicated and not tractable 
because of the presence of the kernel of the Rosenblatt process. 



rl-ff 



m 



oo. 



4 Simulation 

We consider the problem of estimating a from the observations Y±, . . . , Y/v ([!]) driven by the 
approximated Rosenblatt process respectively. In all of the cases, we use a = 2. We have 
implemented the estimator ajy . We have simulated the observations Y\ , . . . , Y N 2 for different 
values of H : 0.55, 0.75 and 0.9 and the values of a: 2 and 20. We consider the cases N = 100 
and N = 200, in others words in the first case we have 10000 observations and in the second 
case 40000. For each case, we calculate 100 estimation djv and we give in the following tables 
the mean and the standard deviation of these estimation. 
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Finally in the model driven by the approximated Rosenblatt process (Section [3D, we 
only construct an estimator that depends on the observations Yj and of the £j . 
The results for N = 100 are the following. 



a = 2 


H = 0.55 


H = 0.75 


H = 0.9 


mean 


2.1860 


2.1324 


2.1416 


stand, deviation 


0.3757 


0.3702 


0.3393 




a = 20 


H = 0.55 


H = 0.75 


H = 0.9 


mean 


20.2643 


20.2013 


20.1506 


stand, deviation 


0.4272 


0.4982 


0.4507 


r N = 200 are the following. 


a = 2 


H = 0.55 


H = 0.75 


H = 0.9 


mean 


2.0433 


2.1048 


2.0559 


stand, deviation 


0.2910 


0.2795 


0.2768 




a = 20 


H = 0.55 


H = 0.75 


H = 0.9 


mean 


20.1895 


20.1080 


20.0943 


stand, deviation 


0.1952 


0.2131 


0.2162 



We can observe that the quality of estimation increases as N increases. We obtain 
these last estimations for the parameter of discretization N = 200. Since the computational 
cost is quite important, we do not implement the estimator for larger N. But although with 
N = 200, the estimator is quite good. 

Appendix: Representation of fBm and Rosenblatt process as 
stochastic integral with respect to a Wiener process 

The fractional Brownian process (-B^)te[o,T] with Hurst parameter H S (0, 1) can be 
written (see e.g. [8], Chapter 5 or [7]) 

Bf = [ K H (t,s)dW s , te[0,T] 
Jo 

where (Wt,t £ [0, T]) is a standard Wiener process, 

K H (t, s) = c H s-2~ H J\u - s) H ~lu H ~^du (12) 

where t > s and ch = [ ^r^wrr i ) 2 an< ^ ^("' ') ^ s ^ ne beta function. For t > s, we have 



dK H . . rs\h-H 



il- 



ls 



An analogous representation for the Rosenblatt process (ZjfOtepVT] i s ( see P2]) 



Zf=f [ F(t, yi ,y 2 )dW yi dW y2 (13) 
Jo Jo 

where (Wt,t £ [0, T]) is a Brownian motion, 

F(t,y 1 ,y 2 ) = d(H)lp it] (yi)l[ 0>t] (y 2 ) / — ^ — — (u,y 2 )du, (14) 



yi Vj/2 5tt 5U 



_ 1 

i?' = — and d(-ff) = jj^i ( 2 (2H-i) ) ■ Actv^rry the representation (fT3ll should be under- 
stood as a multiple integral of order 2 with respect to the Wiener process W. We refer to [8] 
for the construction and the basic properties of multiple stochastic integrals. 
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